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New  energy-based  methods  are  proposed  for  improving  the  performance  of  flexible 
structures  to  external  disturbances.  The  schemes  are  based  on  decentralized  control 
methods,  which  make  them  appealing  for  the  solutions  of  large  structural  control 
problems.  The  proposed  on-off  controller  schemes  use  the  work  done  on  the  actuators 
for  switching.  Energy  consumption  to  damp  out  the  vibrations  is  minimized  by 
switching  the  controllers  ON  only  when  the  work  done  on  them  by  the  structure 
is  positive.  Design  methodology  is  proposed  to  design  the  controllers  to  achieve 
desired  damping  in  vibrating  modes  of  the  structure.  Simulation  results  show  that 
the  proposed  scheme  is  effective  in  introducing  damping  into  low  or  undamped  modes 
of  the  system.  A  method  similar  to  the  Independent  Modal  Space  Control  Scheme  is 
used  to  get  desired  damping  in  modes  of  vibration  using  on-off  control. 


CHAPTER   1 
INTRODUCTION 

1.1     Decentralized  Control  of  Large  Systems 

With  the  space  shuttle  transportation  system  now  a  practical  reality,  there  is 
considerable  interest  in  large  space  structures.  Two  problems  which  are  inherent  in 
the  control  of  large  space  structures  are  shape  control  and  attitude  control.  The  latter 
involves  maintaining  a  given  orientation  of  the  spacecraft  with  regards  to  an  inertial 
reference  system,  while  the  former  involves  the  shape  of  critical  components  of  the 
space  structures  such  as  a  phased  array  antenna.  In  both  cases,  structural  flexibility 
plays  important  roles  since  the  translational  and  attitude  motions  are  coupled  with 
the  structural  vibrations. 

The  motion  of  a  Large  Scale  Flexible  Structure(LFSS)  is  usually  modeled  via 
finite  element  methods,  which  can  result  in  very  large  order  models.  This  means 
that  the  order  of  the  system  may  be  extremely  high,  which  can  make  the  tasks  of 
control  system  design  challenging.  Model  reduction  methods  that  are  applied  in  many- 
cases  (Inman,  1989)  to  control  only  a  subset  of  the  elastic  body  modes  may  lead  to 
spillover  problems,  in  which  the  control  effect  in  stabilizing  the  subset  of  modes  may 
cause  instability  in  the  uncontrolled  modes. 

The  approaches  considered  to  date  for  investigating  the  LFSS  control  have  gen- 
erally been  directed  towards  "centralized  control."  e.g.  model  reduction  methods 
(Hughes  &  Skelton.  1981),  modal  control  methods  (Balas,  1978b),  output  feedback 


control  (Lin.  1981),  and  adaptive  control  methods  (Behnhabib,  I  wens,  k  Jackson, 
1979).  Two  good  survey  papers  deal  with  the  problem  of  LFSS  control  (Gran  & 
Rossi,  1979;  Balas,  1982).  West-vukovich  and  Davison  (1984)  have  shown  how  the 
design  of  active  structural  controllers  that  emulate  the  real  structural  elements  such 
as  dampers  and  springs  can  produce  effective  control  laws.  A  suboptimal  control 
approach  is  used  in  that  work  to  design  the  structural  characteristics  (gains)  of  the 
controller  elements.  This  allows  integration  of  passive  controller,  active  controller 
and  sensor/actuator  locations  design.  This  method  exhibits  problems  due  to  delays 
introduced  due  to  digital  implementation,  and  consideration  of  transfer  function  of 
real  sensors  and  actuators. 

Definition.  Fixed  modes  of  the  system  are  modes  of  the  system  that  remain  in- 
variant under  any  nonzero  parametric  variation  of  the  system  or  the  controller. 

The  fixed  modes  of  a  LFSS  result  when  a  sensor  and  actuator  are  located  at  a 
node  of  a  flexural  mode.  It  was  shown  that  for  LFSS  with  col-located  actuators  and 
sensors,  the  fixed  modes  of  the  decentralized  system  are  the  same  as  the  centralized 
fixed  modes  of  the  system.  A  controller  that  will  eliminate  the  spillover  problem  is 
demonstrated.  It  is  also  shown  that  a  solution  to  the  decentralized  control  problem 
exists  only  if  there  is  a  solution  to  the  centralized  system.  A  further  explanation  of 
fixed  modes  is  provided  in  section  3. 

1.2     Advantage  of  Using  Decentralized  Control 

Over  the  past  twenty-five  years,  engineers  and  scientists  have  developed  a  variety 
of  procedures  for  analyzing  systems  and  for  designing  control  strategies  for  controlling 
LFSS.  These  procedures  are  classified  into  three  types: 

1.  Procedures  for  modeling  dynamical  systems  (state  space  formulation,  input- 
output  transfer  function  description,  etc) 


2.  Procedures  for  describing  the  qualitative  properties  of  system  behavior 
(controllability,  stability,  observability  etc.) 

3.  Procedures  for  controlling  system  behavior  (stabilizing  feedback,  optimal 
control,  etc) 

All  of  these  procedures  rest  on  the  common  presupposition  of  centrality;  all  the  infor- 
mation available  about  the  system,  and  the  calculations  based  upon  this  information, 
are  centralized.  It  is  useful  to  distinguish  between  two  kinds  of  available  information: 

1.  Information  about  the  system  model:  (off-line  or  a  priori  information) 

2.  Sensor  information  about  the  system  response:  set  of  all  real  time  mea- 
surements made  of  system  response. 

The  common  modern  estimation  and  control  is  based  on  centralized  control  the- 
ory. When  considering  large  scale  systems,  the  centralized  control  assumption  fails 
due  either  to  the  lack  of  centralized  information  or  the  lack  of  centralized  comput- 
ing capability.  There  are  many  examples  of  large  scale  systems  that  present  a  great 
challenge  to  both  analysts  and  control  system  designers.  Examples  include  power 
systems,  urban  Traffic  networks,  digital  communication  networks,  flexible  manufac- 
turing networks  and  economic  systems.  The  control  of  such  physical  systems  are  often 
characterized  by  spatial  separation  (such  as  physical  separation  between  sensors  and 
actuators)  so  that  issues  such  as  the  economic  cost  and  reliability  of  the  communi- 
cation link  have  to  be  taken  into  account  in  control,  thus  providing  an  impetus  to 
decentralized  control  design. 

Research  in  decentralized  control  has  been  motivated  by  the  inadequacy  of  modern 
control  theory  to  deal  with  certain  issues  which  are  of  concern  in  large  scale  systems. 
A  key  concept  in  modern  control  theory  is  that  of  state  feedback.  By  using  techniques 
such  as  Linear  Quadratic  (LQ)  optimal  control  or  pole  placement,  it  is  possible  to 


achieve  improved  system  behavior  by  using  state  feedback.  However,  it  is  often  im- 
possible to  instrument  a  system  to  the  extent  required  for  full  state  feedback,  so 
techniques  ranging  from  linear-quadratic  Gaussian  (LQG)  control,  to  observer  based 
control  to  time  domain  compensator  design  techniques  have  been  employed  to  over- 
come this  difficulty.  However,  a  key  concept  in  these  types  of  designs  is  that  every 
sensor  output  affects  every  actuator  input.  This  situation  is  termed  as  centralized 
control.  In  many  large  scale  systems,  it  is  impossible  to  implement  this  in  real-time. 
Thus  for  economic  and  possibly  reliability  reasons,  there  is  a  trend  for  decentral- 
ized decision  making,  distributed  computing  and  hierarchical  control.  However,  these 
desirable  goals  of  structuring  a  distributed  information  and  decision  framework  for 
large  scale  systems  do  not  "mesh"  with  the  available  centralized  methodologies  and 
tools  of  modern  control  system  design.  Reduction  of  computation  and  simplification 
of  structure  are  of  particular  concern  in  decentralized  control  of  large  scale  systems, 
but  are  also  of  concern  in  almost  all  areas  of  control  theory  and  its  applications. 

The  basic  characteristic  of  decentralized  control  is  that  there  are  restrictions  on  the 
information  transfer  between  certain  groups  of  sensors  and  actuators.  For  example, 
in  figure  1,  state  variables  Af  are  used  to  form  the  control  U\  and  state  variables  A'2 
are  used  to  form  the  control  U2-  This  depicts  total  decentralization.  However,  inter- 
mediate restrictions  on  the  information  between  controllers  are  also  possible.  Partial 
decentralization  takes  place  when  the  system  is  not  fully  decentralized  but  the  rate 
of  information  transfer  is  constrained  so  that  full  centralized  control  is  not  possible. 
It  is  to  be  noted  that  the  concept  of  decentralization  refers  to  the  control  structure 
implementation;  the  control  laws  may  be  designed  in  a  completely  centralized  way 
(provided  there  are  no  other  physical  limitations  preventing  this). 

The  advantages  of  using  decentralized  control  were  mentioned  in  the  Ph.D  Thesis 
of  Ahmed  Tarras  and  by  several  other  researchers  (Tarras,  1987;  West-vukovich  & 
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Figure  1.1.  A  decentralized  system 

Davison,  1984).  The  main  advantages  of  using  decentralized  control  laws  are  to 
overcome  the  problems  due  to  the  complexity  of  the  system  that  results  from 

1.  Large  Dimensions 

2.  Uncertainty-  deterministic  or  stochastic 

3.  Structural  Constraints-  these  make  the  flow  of  information  between  the 
subsystems  difficult. 

In  controlling  the  complex  systems,  massive  calculations,  expensive  computer  time 
and  high  computer  costs  may  be  encountered.  In  addition,  there  may  be  large  storage 
requirements,  and  multiple  criteria  may  be  encountered  (due  to  spatial  separations, 
different  levels  of  operation,  etc). 

Definition.  A  control  system  is  called  decentralized,  if  and  only  if  all  local  controls 
are  calculated  only  as  an  explicit  function  of  the  local  information  (states,  output 
etc).  This  is  a  control  with  one  level  and  no  single  controller  has  an  overall  view  of 
the  entire  process. 

In  centralized  control,  one  viable  approach  to  design  of  feedback  control  laws  for 
time  invariant  linear  systems  is  by  minimization  of  an  infinite  horizon  quadratic  per- 
formance index.  LQ  design  method  allows  asymptotic  pole  placement  by  appropriate 
choice  of  the  performance  index  (Kwakarnaak  &  Sivan,  1972)  has  excellent  sensitivity 
and  robustness  properties  according  to  various  criteria,  including  the  classical  gain 


and  phase  margins  (Sofonov  &  Athans.  1977),  and  handles  multi-loop  problems  in 
exactly  the  same  framework  as  single  loop  problems.  These  results  generally  extend 
to  the  stochastic(LQG)  case,  in  which  the  disturbances  are  modeled  by  passing  white 
Gaussian  noise  through  finite  dimensional  filters,  and  in  which  only  partial  state  in- 
formation is  available.  However,  the  solution  of  the  LQG  problem  requires  that  a 
Kalman  filter  be  implemented  to  reconstruct  the  missing  state  variables.  Thus,  an 
additional  Riccati  equation  must  be  solved.  The  modes  of  the  filter  will  appear  in 
the  closed-loop  response  of  the  system,  and  the  gain  and  phase  margins  will  be  more 
sensitive  and  less  robust  than  a  corresponding  LQ  design.  In  nonclassical  stochastic 
control  problems,  where  the  information  available  to  the  device  controlling  a  channel 
input  is  different  from  that  of  the  device  observing  the  channel  output,  the  input 
has  a  dual  purpose:  communication  through  the  system  dynamics  and  sensors  to  the 
other  controllers  and  direct  control  of  the  system. 

From  the  control  sharing  point  of  view,  one  can  speculate  that  nonlinear  solutions 
in  the  nonclassical  LQG  control  are  due  to  actions  of  the  decentralized  controllers  as 
they  try  to  signal  information  to  one  another  using  the  control  system  as  a  commu- 
nication channel.  It  is  impossible  to  ascertain  what  portion  of  a  given  system  input 
is  applied  to  purposes  of  signaling  and  what  portion  for  control.  One  can  question 
the  desirability  of  using  nonlinear  signaling  strategies  to  communicate  through  the 
control  system  dynamics  on  several  grounds.  Implementation  of  such  scheme  would 
be  extremely  complex  and  performance  would  be  quite  sensitive  to  system  parameter 
variations.  In  any  case,  determination  of  these  signaling  strategies  have  been  shown 
to  be  equivalent  to  an  infinite  dimensional,  non-convex  optimal  control  problem  with 
neither  analytical  nor  computational  solution  likely  to  be  forthcoming  (Witsenhausen, 
1973).  Most  research  in  stochastic  decentralized  control  is  concerned  with  determin- 
ing the  optimal  parameters  of  given  decentralized  control  structure  (Date  &  Chow, 
1994;  Levine  &  Athans,  1970). 


1.3     Decentralized  Stabilization  and  Pole  Placement 

A  fundamental  result  in  modern  control  theory  is  that  the  poles  of  a  controllable 
system  can  be  arbitrarily  assigned  (subject  to  complex  pole-pairing  constraints)  by 
state  feedback.  This  result  has  been  extended  to  show  that  the  poles  of  a  closed 
loop  system  consisting  of  controllable  and  observable  linear  system  with  a  dynamic 
compensator  can  be  freely  assigned.  These  results  are  of  great  theoretical  significance 
and  have  served  as  the  basis  of  practical  synthesis  procedure. 

A  natural  generalization  of  the  pole-placement  question  arises  when  the  restriction 
to  decentralized  feedback  control  is  made.  Although  several  authors  had  looked  at 
this  question,(McFadden,  1967;  Aoki,  1972;  Corfmat  &  Morse.  1973)  the  most  distinct 
results  are  those  of  Wang  and  Davison  (1973),  Davison  (1976)  and  R.P.Corfmat  and 
A.S.Morse  (1976).  Here,  their  results  are  briefly  summarized. 

For  a  linear  system  the  problem  of  decentralized  pole  placement  can  be  formulated 
as  follows.  Consider  the  linear  system 


jv 
x{t)  =  Ax{t)+YdBiui{t)  (1.3.1) 


y,(t)  =  C,x(t)  (1-3-2) 

where  i  =  1 N  indexes  the  input  and  output  variables  of  the  various  controllers, 

x  €  K"  is  the  state,  u;  e  Rm>  and  </,  G  RT<  are  the  input  and  output,  respectively,  of 
the  i'h  local  control  station.  .4  is  the  state  matrix,  Bi  is  the  input  matrix  and  C,  is  the 
output  matrix  of  control  station  i.  The  itk  controller  employs  dynamic  compensation 
of  the  form 

Ui(t)  =  MiZi(t)  +  Fiyt(t)  +  GMt)  (1.3.3) 

~,(t)  =  HiZiit)  +  LlVl(t)  +  RiVt(t)  (1.3.4) 


where  z,  6  Rm  is  the  state  of  the  ith  feedback  controller,  vt  6  Rm<  is  the  ith 
local  external  input,  and  M,.  Ff,  Gu  Hit  L{,  Rt  are  real  constant  matrices  of 
appropriate  sizes.  The  decentralized  pole  placement  problem  is  to  find  matrices 
Mi,Fi,Gi.Hi,L,.Rt  such  that  the  closed  loop  system  described  by  (1.3.1)-(1.3.4) 
has  prespecified  poles.  Of  course,  if  [C„  A,  Bt\  is  controllable  and  observable  from 
all  the  stations,  the  solution  to  such  problem  is  guaranteed.  The  interesting  case  is 

to  assume  that  (1.3.2)  is  controllable  from  all  controls  ut uN  but  not  from  any 

single  control  u,,  with  similar  observability  assumptions. 

Consider  first  the  special  case  Mt  =  0  in  the  above  problem.  This  corresponds 
to  a  static  decentralized  output  feedback  controller.    If  F  denotes  the  collection  of 

feedback  matrices  (Fi,F2, Fy),  then  the  pole  placement  problem  is  to  determine 

F  such  that  the  matrix 

Af^A  +  ^BtFtQ  (1.3.5) 


has  an  arbitrarily  specified  set  of  eigenvalues.  Clearly,  a  necessary  condition  for  pole 
placement  in  this  case  is  that  the  polynomials  [  XI  -  AF  |  have  no  common  factor, 
i.e.  that 

a(X)  =  g.c.d.  \XI-AF\=l  (1.3.6) 

where  g.c.d.  is  greatest  common  divisor.  What  is  more  interesting  is  that  this  con- 
dition is  both  necessary  and  sufficient  for  pole  placement  with  dynamic  compensa- 
tion(Wang  &  Davison,  1973).  More  generally,  since  the  zeros  of  a  (A)  (termed  as  the 
fixed  modes  of  the  system)  are  invariant  under  decentralized  dynamic  compensation, 
it  follows  that  a  necessary  and  sufficient  condition  for  stabilizability  is  that  the  roots 
of  a  (A)  have  strictly  negative  real  parts. 

Computation  of  the  fixed  modes  of  a  system  can  be  accomplished  by  computing 
the  eigenvalues  of  AF  and  .4  for  randomly  selected  F  and  checking  for  common 
eigenvalues.  A  simpler  way  of  checking  for  the  fixed  modes  of  the  system  is  given  in 


Anderson  and  Moore  (1981).  Implicit  in  the  pole  placement  result  quoted  above  is 
a  constructive  algorithm.  This  algorithm  requires  as  a  first  step  the  selection  of  F 
such  that  the  poles  of  AF  are  distinct  from  that  of  A.  Then,  the  dynamic  feedback  is 
successively  employed  at  the  control  stations  to  place  the  poles  that  are  controllable 
and  observable  from  a  given  station. 

R.P.Corfmat  and  A.S.Morse  (1976)  have  studied  the  decentralized  feedback  con- 
trol problem  from  the  point  of  view  of  determining  a  more  complete  characterization 
of  conditions  for  stabilizablity  and  pole  placement.  Their  basic  approach  is  to  de- 
termine conditions  under  which  a  system  of  the  form  (1.3.1)-(1.3.4)  can  be  made 
controllable  and  observable  from  the  input  and  output  variables  of  a  given  system 
by  applying  static  feedback  to  other  controllers.  Then  dynamic  compensation  can  be 
employed  at  this  controller  in  a  standard  way  to  place  the  poles  of  the  system. 

It  is  not  hard  to  see  that  a  necessary  and  sufficient  condition  to  make  (1.3.1)- 
(1.3.4)  controllable  and  observable  from  a  single  controller  is  that  none  of  the  transfer 
functions 

Gy{8)  =  Ci(sl  -  A)-lBj  i,j  =  l, ,N  (1.3.7) 

vanish  identically.  A  system  satisfying  this  condition  is  termed  strongly  connected. 

If  a  system  is  not  strongly  connected,  it  is  impossible  to  make  the  system  con- 
trollable and  observable  from  a  single  controller.  In  this  case,  it  is  necessary  to 
decompose  the  system  into  a  set  of  strongly  connected  subsystems  and  to  make  each 
subsystem  controllable  and  observable  from  one  of  its  controllers.  For  a  strongly 
connected  system,  Corfmat  and  Morse  have  given  a  highly  interesting  and  rather  in- 
tuitive condition  that  is  necessary  and  sufficient  to  make  (1.3.1)-(1.3.4)  controllable 
and  observable  from  a  single  controller.  They  have  shown  that  if  a  strongly  connected 
system  can  be  made  controllable  and  observable  from  a  single  controller,  it  can  be 
made  controllable  and  observable  from  anv  controller,  and  a  necessary  and  sufficient 
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condition  for  this  to  occur  is  that  the  system  be  complete.  Completeness  of  the  sys- 
tem is  defined  in  terms  of  the  transmission  polynomials  of  the  system.  For  example, 
in  case  N  =  2,  a  test  for  completeness  (Morse,  1973)  is  that  for  all  s 


rank 


(si -A)    B, 
Co  0 


>  2,  (1.3.8) 


rank  (si  -  A)  >  2  (1.3.9) 

except  when  s  is  an  eigenvalue  of  A.  It  is  only  necessary  to  check  the  above  condition 
on  the  spectrum  of  A. 

Anderson  and  Moore  (1981)  suggested  that  time  varying  controllers  can  be  used  to 
eliminate  the  fixed  modes.  In  their  work,  they  showed  that  under  the  assumption  of 
centralized  controllability  and  observability,  for  a  two  channel  system  with  feedback 
of  the  form:  u2(t)  —  K2(t)y2(t)  where  K2(t)  is  periodic  and  piecewise  constant,  taking 
p  >  1  +  max  (dim(u2) ,  dim(y2))  values,  then  strong  connectivity,  even  when  a  fixed 
mode  is  present,  is  enough  to  ensure  that  the  system  is  uniformly  controllable  from 
u\  and  uniformly  observable  from  yi. 

Most  of  the  work  in  the  decentralized  control  schemes  with  time  varying  controller 
focuses  on  making  the  system  controllable  and  observable  and/or  on  robust  control. 
The  concept  is  to  develop  an  algorithm  to  place  the  poles  of  the  systems  using  time 
varying  controllers.  Once  this  algorithm  is  available,  fixed  modes  (  the  modes  of  the 
system  which  are  not  affected  by  TIV  controller  gains)  of  the  system  can  be  treated 
the  way  the  non-fixed  modes  are  treated  (Wang,  1982). 

The  concept  of  using  time  varying  periodic  compensators  with  time  invariant 
plants  was  further  investigated  by  Khargonekar,  Poolla,  and  Tannenbaum  (1985).  In 
their  work,  they  showed  that  for  a  large  class  of  robustness  problems,  periodic  com- 
pensators are  superior  to  time  invariant  ones.  They  also  gave  design  techniques  which 
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can  be  implemented  easily.  They  also  showed  that  for  weighted  sensitivity  minimiza- 
tion for  linear  time  invariant  plants,  time  varying  controllers  have  no  advantages  over 
time  invariant  ones. 

However,  the  eigenstructure  assignment  problem  was  beyond  the  scope  of  Poolla's 
work.  Their  design  was  used  to  design  a  stabilizing  compensator  and  later,  the  time 
varying  controller  and  time  invariant  plant  combinations  were  used  to  generate  the 
time  response  of  the  system.  Identification  packages  were  used  to  identify  this  system 
and  separate  compensator  was  designed  for  pole  placement.  The  order  of  the  system 
seems  to  be  going  higher  and  higher  in  this  case.  Nevertheless,  this  looks  like  a  first 
pass  solution  to  the  problem  of  eigenstructure  assignment. 

1.4     An  illustration  to  demonstrate  the  advantages  of  decentralized  structural  control 

This  example  shows  the  advantages  and  simplicity  of  building  a  decentralized 
controller  for  structures.  Figure  1.4  shows  a  6-bay  truss  whose  model  is  taken  from 
Young  (1990). 

This  truss  is  modeled  as  a  coupled  structure  of  three  substructures.  Substructure 
boundaries  are  marked  in  the  figure.  The  truss  member  mass  and  stiffness  matrices 
expressed  with  respect  to  the  local  coordinates  and  used  in  the  assembly  process  are 
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6 


2     1 

1    2 


The  nodal  coordinates  are  defined  as  vertical  and  horizontal  displacements  at  the 
joints.  The  internal  degrees  of  freedom  at  which  the  actuator  and  displacement  sensors 
are  placed  are  marked  in  the  figure.  The  finite  element  model  of  the  individual  sub- 
structures are  obtained  using  a  FORTRAN  program  and  validated  using  the  common 
FEM  package  ANSYS  (Swanson  Analysis  Inc,  1992).  For  convenience,  the  material 
properties  are  assumed  to  be  have  unit  magnitudes.   An  Interlocking  control  design 
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concept  is  used  to  place  the  actuators  and  sensors  at  internal  degrees  of  freedom.  Min- 
imization of  internal  coordinate  motions  of  different  substructures  would  localize  the 
dynamic  interaction  of  coupled  structure  in  the  components.  The  component  control 
action  is  designed  to  lock  up  its  own  boundary  condition,  that  better  approximates 
the  one  assumed  in  the  component  modeling  of  its  adjacent  components. 


Figure  1.2.  A  Six-Bay  Truss 

A  convenient  control  design  technique  for  this  concept  is  the  linear  quadratic  op- 
timal control  regulator  approach,  in  which  the  aim  is  to  minimize  the  internal  coordi- 
nates motion.  The  component  control  law  is  the  one  that  minimizes  the  performance 
index. 


Ji  =  \  /o°°  {ysTys  +  wTR*f)  dt 


where  ys  and  u'  denote  the  output  and  control  inputs  at  the  internal  degrees  of 
freedom  of  each  substructure.  These  are  marked  in  the  figure  by  numbers.  The  opti- 
mal component  control  is  a  state  feedback  control  law  with  positions  and  velocities  as 
inputs.  This  type  of  controller  is  built  for  individual  substructures  and  the  controlled 
coupled  structures  state  matrix  is  obtained.  The  effect  of  having  such  a  controller 
on  the  eigenvalues  is  studied.  Failures  in  the  component  controllers  by  changing  the 
elements  of  the  gain  matrices  of  individual  substructures.  The  gains  of  some  channels 
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were  set  to  zero  (  a  scenario  mimicking  a  break  in  communication  line)  and  the  effect 
of  this  on  the  overall  eigenvalues  were  studied. 

With  the  decentralized  control,  it  was  observed  that  all  the  poles  of  the  coupled 
structure  have  negative  real  parts.  The  individual  controller  failures  don't  seem  to 
affect  the  pole  locations  very  much.  This  of  course  depends  on  the  type  of  controller 
failure  that  we  simulate.  Open  loop  poles  of  the  substructures  and  coupled  structure 
are  plotted  in  Figure(1.4).  In  table  1.  Values  of  the  coupled  loop  poles  without  and 
with  controller  failures  are  listed. 
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Figure  1.3.  Open  loop  poles  of  substructures  and  coupled  structure 


1.5     Objective  of  the  Present  Study 


The  present  study  investigates  the  development  of  new  and  promising  energy  effi- 
cient ways  of  controlling  vibrations  in  mechanical  structures.  Although  considerable 
research  has  been  done  in  vibration  control  of  flexible  structures,  no  method  uses 
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Without  Component 

With  End  controlled 

With  Center  Controlled 

failure 

Component  Failure 

Component  Failure 

-0.0005  +  0.1877i 

-0.3e-8  +  0.1876i 

-0.0005  +  0.18771 

-0.0005  +  0.1877i 

-0.0005  +  0.1877i 

-0.0005  +  0.1877i 

-0.0037  +  0.2515i 

-0.0037  +  0.2515i 

-0.0003  +  0.2514i 

-0.0007  +  0.2679i 

-0.0004  +  0.2679i 

-0.0003  +  0.2679i 

-0.0101  +  0.6109i 

-0.0004  +  0.6092i 

-0.0097  +  0.61061 

-0.0093  +  0.61171 

-0.0093  +  0.6117i 

-0.0093  +  0.6117i 

-0.0054  +  0.6437i 

-0.0024  +  0.6428i 

-0.0035  +  0.644H 

-0.0059  +  0.64451 

-0.0042  +  0.6448i 

-0.0027  +  0.644H 

-0.0044  +  0.6488i 

-0.0036  +  0.6489i 

-0.0030  +  0.6485i 

-0.0029  +  0.6563i 

-0.0028  +  0.6563i 

-0.0023  +  0.6562i 

-0.0059  +  0.6909i 

-0.0001  +  0.6903i 

-0.0058  +  0.6907i 

-0.0055  +  0.6910i 

-0.0056  +  0.6910i 

-0.0055  +  0.6910i 

-0.0196  +  1.03431 

-0.0181  +  1 .0343i 

-0.0041  +  1.03091 

-0.0062  +  1.0600i 

-0.0029  +  1.05871 

-0.0047  +  1.0596i 

-0.0045  +  1.2389i 

0.0000  +  1.2375i 

-0.0044  +  1.2388i 

-0.0043  +  1.2389i 

-0.0044  +  1.2389i 

-0.0043  +  1.2389i 

-0.6676  +  2.32061 

-0.6674  +  2.3318i 

-0.6991  +  2.4474i 

-0.6813  +  2.39331 

-0.6932  +  2.4306J 

-0.6997  +  2.45571 

-0.7125  +  2.4606i 

-0.7125  +  2.4642i 

-0.6914  +  2.4778i 

-0.7252  +  2.4778i 

-0.0058  +  2.54181 

-0.0126  +  2.48181 

-0.7066  +  2.4807i 

-0.7247  +  2.5627i 

-0.7084  +  2.482H 

-0.7325  +  2.5748i 

-0.0044  +  2.5898i 

-0.0162  +  2.5758i 

-1.1823  +  2.8937i 

-1.1842  +  2.9214i 

-1.2102  +  2.9729i 

-1.2108  +  2.9733i 

-1.2149  +  2.978H 

-1.2061  +  2.97791 

-1.2238  +  3.0177i 

-0.0079  +  3.2286i 

-0.0085  +  3.1672i 

-2.2575  +  3.76781 

-2.2922  +  3.7970i 

-2.3014  +  3.8008i 

-2.3022  +  3.801  li 

-2.3427  +  3.859H 

-2.3012  +  3.80401 

-2.3828  +  3.88871 

-0.0032  +  4.4488i 

-0.0067  +  4.5129i 

Table  1.1.  Eignevalues  of  the  6-bay  truss  example 


energy  transfer  from  the  structure  to  the  controller  as  a  criterion  to  decide  on-off 
times.  The  main  objective  is  to  demonstrate  that  such  a  method  is  viable  and  to  give 
a  controller  design  methodology. 

In  Chapter  2.  analytical  results  are  presented  to  demonstrate  the  use  of  on-off 
controllers  to  damp  the  vibrations  in  an  energy  efficient  way.  Energy  input  into  the 
controller  is  used  as  a  criterion  to  select  the  switching  times.  The  controllers  are 
switched  on  only  when  work  is  done  on  them.    This  does  not  require  the  use  of  an 
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external  energy  source.  Since  energy  supplied  to  the  controller  would  require  that 
fuel  be  consumed,  the  proposed  design  methodology  reduces  the  fuel  consumption  for 
vibration  suppression.  Additionally,  the  proposed  control  laws  are  decentralized  since 
the  on-off  switching  laws  are  based  on  the  local  displacements  and  velocities. 

In  Chapter  3.  controller  design  methodologies  are  given  to  select  the  controller 
parameters  to  achieve  the  desired  modal  damping.  Analytical  results  are  given  to 
prove  the  validity  of  this.  A  method  which  is  similar  to  the  Modified  Independent 
Modal  Space  Control  (MIMSC)  is  used  for  designing  on-off  controller  for  multi  degree 
of  freedom  systems.  Simulation  results  are  presented  to  show  the  effectiveness  and 
accuracy  of  this  method. 

In  Chapter  4,  the  use  of  optimal  control  laws  with  fuel  consumption  as  a  cost 
criterion  to  design  vibration  controllers  is  presented.  First,  the  concept  of  on-off 
vibration  control  is  reviewed  with  an  emphasis  on  fuel  consumption  minimization. 
Next,  the  problems  associated  with  window  sizing  and  selecting  the  sampling  interval 
are  brought  to  light. 

A  discussion  of  energy  efficient  switching  control  laws  future  research  problems 
are  presented  in  Chapter  5. 


CHAPTER   2 
MINIMUM  ENERGY  ON-OFF  CONTROL  FOR  MECHANICAL  SYSTEMS 

2.1     Introduction 

In  this  work  a  minimum  energy  control  method  for  on-off  decentralized  control 
of  mechanical  systems  is  introduced.  Energy  consumption  is  minimized  by  turning 
on  the  controllers  when  the  structure  does  work  on  the  controller  and  turning  them 
off  when  the  controller  would  impart  energy  on  the  structure.  Vibrating  energy  from 
lightly  damped  modes  is  transferred  to  highly  damped  modes  to  introduce  damping. 
In  the  next  chapter,  this  method  is  applied  to  solve  pole  placement  problems  for 
vibrating  systems. 

The  work  presented  here  is  based  on  the  intuition  that  using  an  energy  minimiza- 
tion criterion  that  works  with  conservation  principles  will  produce  energy  efficient 
methods  for  control  of  vibrating  systems.  The  methodology  presented  here  is  quite 
feasible,  since  the  switching  of  the  electro-mechanical  controller  is  performed  electron- 
ically via  an  on-board  processor.  The  development  uses  energy  as  an  optimization 
criterion.  The  delays  introduced  by  the  decision  making  are  not  included,  but  can  be 
very  easily  incorporated  into  the  analytical  results  presented  here. 

Several  authors  used  or  proposed  the  use  of  on-off  control  for  attitude  and  shape 
control  of  large  flexible  space  structures  (  Yelde  and  He  (1983),  Foster  and  Silver- 
berg  (1991),  Arbel  and  Gupta  (1981).  Masri.  Bekey,  and  Caughey  (1981),  Rohman 
and  Leipholz  (1978),  Seywald,  Kumar.  Deshpande.  and  Heck  (1994).  Redmond  and 
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Silverberg  (1992).  Neustadt  (1960)  and  Reinhorn.  Soong,  Riley,  Lin,  Aizawa.  and 
Higashino  (1993)).  The  approaches  using  on-off  control  devices  are  relatively  simple 
and  require  less  on-line  computational  effort  than  other  modern  control  techniques. 
High  energy  pulses  can  be  input  into  the  system  using  on-off  devices.  In  the  works 
of  Zimmerman.  Inman,  and  Juang  (1991),  Zimmerman  (1990),  Masri  et  al.  (1981),  in 
order  to  conserve  energy,  control  devices  are  activated  when  some  specified  thresholds 
of  states  has  been  exceeded.  The  pulse  magnitude  was  determined  analytically  so  as 
to  minimize  a  non-negative  cost  function  related  to  the  system  energy.  The  degree 
of  control  obtained  depended  on  the  threshold  level  considered  and  the  cost  function 
that  was  minimized. 

Foster  and  Silverberg  (1991),  Arbel  and  Gupta  (1981),  Masri  et  al.  (1981),  Sey- 
wald  et  al.  (1994),  Redmond  and  Silverberg  (1992)  and  Neustadt  (1960)  used  the 
principle  of  fuel  minimization  as  a  criterion  for  the  implementation  of  on-off  switch- 
ing control.  Flight  duration  of  spacecraft  is  limited  in  many  cases  by  the  amount 
of  on-board  fuel.  For  satellites,  on-board  fuel  is  needed  to  fire  the  thrusters  and 
apogee/perigee  motors  to  place  them  in  the  final  orbit.  Thruster  firings  will  be  needed 
for  station-keeping  and  to  correct  orbital  errors.  As  the  size  of  spacecraft  increases, 
the  need  to  suppress  the  vibrations  introduced  into  the  spacecraft  body  must  also  be 
incorporated  into  the  controller  design. 

The  work  described  in  this  chapter  differs  from  previous  works  in  that,  here,  fuel 
consumption  is  measured  by  the  amount  of  external  energy  needed  to  damp  out  the 
vibrations.  More  discussions  about  the  proposed  scheme  and  the  types  of  controllers 
that  may  be  used  with  the  proposed  scheme  may  be  found  in  the  next  two  sections. 

2.2     Controlling  Vibrations  Using  On-Off  Controllers  Without  External  Energy 

In  general,  when  subject  to  vibrational  control,  a  structure  does  work  on  its  con- 
troller and  vice- versa.  Here  it  is  proposed  that  the  work  done  on  the  controllers  by  the 
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vibrating  structure  be  taken  advantage  of  to  arrive  at  a  fuel-optimal  solution.  This  is 
accomplished  by  switching  the  controller  "on"  only  when  the  structure  is  doing  work 
on  the  controller  and  switching  it  ''off"  otherwise. 

2.2.1     Switching  procedure 

In  this  section  it  is  demonstrated  how  an  on-off  control  scheme  using  only  the 
energy  input  rate  into  the  controller  can  get  a  better  response  from  the  system. 
Preliminary  versions  of  this  work  are  presented  in  Schwartz  and  Maben  (1996.  1995). 

Consider  an  n-dimensional  second  order  dynamical  system  of  the  form 

Mx  +  Dx  +  Kx  =  0.  (2.2.1) 

where  M,  D,  K  are  the  mass,  damping  and  stiffness  matrices  respectively.  D  is 
assumed  to  be  of  proportional  type,  i.e.  D  satisfies  D  =  aM  +  0K,  where  a,  0  are 
scalar  constants.  Also,  K  =  KT  >  0  and  M  =  MT  >  0. 

When  the  controller  is  off,  the  energy  in  the  system  is  given  by 

EoU  =  -{xTMx  +  xTKx  -  xTDx\.  (2.2.2) 

Since  the  system  is  dissipative,  the  rate  of  change  of  energy  is  given  by 

EBf,  =  xT[Mx  +  (A»  -  xTDx.  (2.2.3) 

If  the  controller  is  of  the  spring  type,  the  rate  of  change  of  energy  while  the 
controller  is  on  is  given  by 

Em  =  xT[Mx  +  (K  +  AA'ji]  -  xTDx.  (2.2.4) 

where  AA  is  the  added  stiffness  due  to  the  controller.  It  is  assumed  that  AK  is 
symmetric  and  positive  definite. 

In  this  study  we  are  not  concerned  with  measuring  the  energy  stored  in  a  con- 
trolling device.    Rather  we  are  interested  in  studying  the  energy  input  into  such  a 
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device  by  observing  the  rate  of  change  of  energy  function.  Our  aim  is  to  increase  the 
energy  into  such  devices  (on  the  assumption  that  a  unidirectional  force  on  a  spring 
type  of  controller  can  be  converted  into  some  other  form  of  energy  and  released  into 
an  energy  dissipating  or  storage  device  and  not  released  back  into  the  structure). 
The  controllers  are  thus  switched  on  only  when  this  energy  rate  is  positive.  It  should 
also  be  noted  that  x  in  all  these  cases  represents  the  displacement  over  a  bias  value, 
measured  in  reference  to  an  inertial  level. 

The  signal  Eon  -  E„/f  =  xTAKx  gives  a  measure  of  the  rate  of  work  done  on 
all  controllers  if  they  were  on  continuously.  In  the  case  of  one  controller,  if  AA'„  is 
the  stiffness  introduced  by  the  controller,  xa  is  the  relative  displacement  and  xa  is 
the  velocity  at  the  controller  location,  Em  —  E0ff  is  positive  (work  is  done  on  the 
controller  by  the  structure),  if  the  product  of  force  acting  on  the  controller  AKaxa 
and  the  velocity  at  the  controller  location  xa  is  positive.  In  the  case  of  multiple 
controllers,  if  i  denotes  the  location  of  the  itk  controller,  [AAji;]  is  monitored  and 
each  controller  is  switched  on  when  the  force  signal  acting  on  it  and  the  velocity 
at  the  controller  location  are  positive.  The  force  acting  on  the  elements  of  AK 
matrix  can  be  positive  or  negative.  Assuming  that  the  force  in  the  positive  direction 
(representing  compression  in  spring  type  of  controller)  can  be  released  into  some  other 
energy  dissipating  device  (for  example  a  resistor  network),  controllers  are  switched 
on  only  when  the  forces  acting  on  them  are  positive.  This  concept  is  summarized  in 
the  following  switching  procedure. 

It  should  be  noted  that  the  location  of  the  controller  plays  a  role  in  the  controlla- 
bility of  the  vibrating  modes  of  the  structure.  If  the  location  of  an  controller  is  close 
to  the  node  of  a  vibrating  mode,  that  controller  will  have  very  little  effect  on  the 
damping  of  that  particular  mode  of  vibration. 
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2.2.2  Procedure  for  switching  with  one  controller 

1.  Between  switching  times,  monitor  the  signals  AKaxa  and  xa  where  AKa  is  the 
nonzero  entry  on  the  AA'  matrix. 

2.  When  both  AKaxa  and  xu  are  positive,  switch  the  controller  on.  Switch  the 
controller  off  at  the  end  of  a  predetermined  time  which  is  related  to  the  system 
time  constant  as  determined  by  the  lowest  period  of  oscillation  in  the  structure.1 
(and  release  the  energy  absorbed  by  the  spring  when  it  is  switched  off). 

2.2.3  Example 

Consider  a  5  degree  of  freedom  spring  mass  system  shown  in  Figure  2.2.3.  Here  we 
have  a  second  order  system,  Mx+Cx+Kx  =  /  ,  where  the  mass  and  stiffness  matrices 


are  chosen  as,  mass  = 
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Figure  2.1.  A  5  DOF  Spring/Mass  System 


In  this  example  xa  is  X4.  All  but  the  fourth  mode  have  very  little  damping  (equiv- 
alent damping  factor  of  0.002),  while  the  fourth  mode  is  highly  damped(equivalent 
damping  factor  of  0.4).  The  simulations  of  the  control  procedure  of  Section  2.2.2 
show  this  method  to  be  very  effective  in  damping  out  the  vibrations  of  the  measured 


LIn  the  example  in  the  next  subsection,  this  time  is  set  to  be  one  tenth  of  lowest  period  of 
oscillation. 
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response  even  when  the  initial  conditions  were  very  close  to  the  mode  shapes  of  a 
lesser  damped  mode  of  the  uncontrolled  structure,  a  worst  case  scenario. 

Figures  2.2  through  2.6  show  the  nodal  (position)  response  of  the  uncontrolled 
and  controlled  system  when  the  initial  conditions  were  set  to  the  mode  shape  of  one 
of  the  lesser  damped  mode  of  the  uncontrolled  structure.  The  results  are  only  shown 
for  one  set  of  initial  conditions,  however  the  simulation  was  carried  out  for  each  of 
the  mode  shapes.  From  these  plots  it  is  clear  that  the  proposed  method  is  effective 
in  increasing  the  damping. 


Figure  2.2.  Response  at  node  1 
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Figure  2.3.  Response  at  node  2 


Figure  2.4.  Response  at  node  3 
2.3     Summary 


A  minimum  energy  on-off  vibration  control  law  for  structures  was  introduced. 
The  control  procedure  takes  advantage  of  the  energy  which  can  be  imparted  to  the 
controllers  from  switching  by  turning  on  the  controllers  when  the  structure  is  doing 
work  on  the  controller  and  turning  them  off  when  the  controllers  are  doing  work  on 
the  structure. 


23 


iliiiiliiliiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiil'i 


Figure  2.5.  Response  at  node  4 


Time  in  seconds 


Figure  2.6.  Response  at  node  5 


The  proposed  scheme  uses  the  fact  that  a  structure  does  work  on  the  controllers 
and  vice  versa  while  undergoing  vibrations.  The  controllers  are  switched  on  only 
when  they  are  absorbing  energy  from  the  structure.  This  absorbed  energy  can  be 
transformed  into  some  other  form  (to  do  some  useful  work  like  storing  energy  in  a 
device  like  battery  or  gyro).  The  proposed  scheme  achieves  the  desired  goal  without 
the  use  of  an  external  power  source  (except  for  decision  making  and  measurements). 
In  one  application  which  gives  promise  to  the  proposed  scheme,  Yang,  Mikulas.  Park, 
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and  Su  (1995)  have  shown  that  controlled  moment  gyros  can  be  used  for  the  slewing 
maneuvers  and  vibration  control  of  space  structures. 


CHAPTER  3 
DESIGNING  CONTROLLER  TO  ACHIEVE  DESIRED  DAMPING 

3.1     Introduction 

This  chapter  provides  methods  for  choosing  controller  parameters  AK  in  order  to 
achieve  specified  modal  damping  for  single  degree  of  freedom  systems  (SDOF)  and 
multi  degree  of  freedom  systems  (MDOF).  First,  analytical  results  are  presented  for 
SDOF  systems.  Second,  results  are  given  for  MDOF  systems.  Modified  Independent 
Modal  Control  Methods  (MIMC)  are  applied  to  MDOF  systems  to  achieve  the  desired 
modal  damping.  It  is  also  shown  that  introducing  damping  in  only  few  modes  of  the 
uncontrolled  system  may  cause  poor  responses  in  other  modes.  In  view  of  this,  it  is 
recommended  that  all  dominant  modes  of  the  structure  be  damped  when  using  the 
proposed  switching  control  design. 

3.2     Designing  On-Off  Controller  for  SDOF  system 

3.2.1     Analytical  Results 

Consider  the  scalar  single  of  freedom  system  described  by 

mx  +  kx  +  fiAkx  =  0,  (3.2.1) 

where  m  is  the  mass,  A:  is  the  stiffness,  A/c  is  the  added  stiffness  due  to  the  controller 
action,  x  is  the  displacement  of  the  mass,  and  /i  is  a  parameter  which  takes  values  1 
or  0,  depending,  respectively,  upon  whether  or  not  the  controller  is  on.  The  switching 
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control  is  expected  to  introduce  damping  into  the  system.  The  controller  is  assumed 
to  be  located  at  the  node  associated  with  this  mass  and  Ak  is  positive.  The  control 
objective  is  to  achieve  a  specified  average  damping  rate.   The  following  is  a  design 
procedure  for  specifying  Ak  so  as  to  achieve  this  objective. 
The  energy  in  the  system  is  given  by 

E  =  \  [mx2  +  kx2  +  ftAkx2)  .  (3.2.2) 

Let  w  be  the  natural  frequency  of  the  uncontrolled  system,  (i.e.  w2  =  £j  . 
The  transformation  \/mx  =  u  gives 

E  =  ~  U2  +  wV  +  /^2=^'«2)  •  (3-2-3) 

By  imposing  the  control  scheme  described  in  the  previous  chapter,  fi  is  set  to  1 
when  t  is  such  that  Aku  >  0,  or  equivalently  when  u  >  0. 

It  is  expected  that  the  proposed  switching  control  scheme  introduces  damping  so 
that  u  is  assumed  to  be  of  the  form: 

u  =  asm{ut  +  p),  (3.2.4) 

where  a  =  a(t)  is  an  exponentially  decaying  function  of  time,  and  w  is  the  uncon- 
trolled natural  frequency  of  oscillation.  A  functional  relationship  between  the  rate 
of  decay  function  a  and  the  feedback  gain  Ak  will  be  used  to  determine  a  law  for 
introducing  a  desired  average  damping  into  the  system. 

Taking  the  time  derivative  of  both  sides  of  equation  (3.2.4), 

ii  =  auJCOs(u)t  +  0)  +  cesin(ujt  +  0).  (3.2.5) 

The  decay  rate  a  is  assumed  to  be  small  compared  to  aw,  so  the  a  term  may  be 
neglected  in  equation  (3.2.5).  This  gives 
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ii  zx  auj  cos{wt  +  rf) .  (3.2.6) 

By  conservation  of  energy,  E  =  0.  Note  that  this  is  true  at  all  instances  except 
switching:  Thus 

1  d  Ak   \ 
— -(it2  +  uj'2u2)  +  <hiuj2—u     =  0.  (3.2.7) 

2  at  k    J 

Using  equations  (3.2.6)  and  (3.2.4)  in  equation  (3.2.7)  yields 

r,  UU2     Ak  1 

i02aa  +  ^- — —aiusm2(ut  +  P)  =  Q.  (3.2.8) 

2      A: 

Since  the  decay  rate  of  the  response  is  assumed  to  be  slow,  the  method  of  averaging 

in  (Caughey  &  O'Kelly,  1965)  may  be  applied  to  compute  d: 

Ak    uj    r11/^ 
q= — Ta7;  iism2(ut  +  3)dt.  (3.2.9) 

k      2  Jo 

Since  u  >  0  holds  for  only  half  a  cycle  of  sin2(u;£  +  8),  fi  is  set  to  1  for  only  half 
of  a  cycle  of  sin2(wi  +  0).  In  that  case,  the  value  of  the  integral  in  equation  (3.2.9) 
is  1/ui  and 


(3.2.10) 


Ak 
a  =  Q0e_<:t.     where    £  =  — .  (3.2.11) 

Thus,  for  the  single  degree  of  freedom  case,  determines  a  control  law  for  achieving  a 
desired  average  damping  rate  £  =  $-. 

This  result  is  summarized  in  the  following  proposition. 


Ak      w 

Ak 

d  = 

k      2ttui 

k 
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3.2.2  Proposition 

For  a  single  degree  of  freedom  system  of  the  form  (3.2.1),  on-off  switching  control 
may  be  used  to  attain  a  desired  overall  average  damping  rate  (  using  the  proposed 
minimum  energy  control  scheme  in  conjunction  with  equation  (3.2.11) 

3.2.3  SDOF  Example 

The  main  objective  of  this  example  is  to  show  the  effectiveness  of  this  method  and 
to  show  the  validity  of  assumption  made  earlier,  that  the  decay  rate  a  is  assumed  to 
be  small  compared  to  auj. 

Consider  the  single  degree  of  freedom  system  with  magnitudes  of  mass  and  stiff- 
ness respectively  given  by  m=l  and  k=100.  Simulation  of  the  system  with  on-off 
controller  was  carried  out  for  three  values  of  the  damping.  For  each  of  these  damp- 
ing, Ak  was  calculated  using  equation  (3.2.11).  The  on-off  controller  was  simulated 
and  the  response  of  the  system  to  an  initial  displacement  of  1  and  velocity  of  0  was 
observed.  The  response  of  the  system  for  values  of  damping  of  0.0015,  0.020  and 
.6  are  given  in  Figures  3.1,  3.2  and  3.3  respectively.  In  each  of  these  figures,  expo- 
nentially decayed  curves  are  plotted  to  show  the  exactness  of  this  method  for  low 
values  of  damping.  From  these  time  responses,  it  is  concluded  that  the  assumption 
that  a  is  an  exponentially  decayed  function  is  valid.  Also,  it  can  be  seen  that  when 
the  assumption  that  a  «  aw  is  violated,  the  damping  introduced  into  the  system 
because  of  switching  controller  deviates  much  from  the  desired  value. 
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Figure  3.1.  Response  of  SDOF  £  =  0.0015 
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Figure  3.2.  Response  of  SDOF  £  =  0.02 
3.3     Designing  On-Off  Controll  for  an  MDOF  system 
3.3.1     Modified  Independent  Modal  Control  Methods 


Active  control  of  the  vibrations  of  flexible  structures  based  primarily  on  modal 
control  methods  whereby  the  vibrations  are  suppressed  by  controlling  the  dominant 
modes  of  vibration  has  been  the  focus  of  many  researchers  (Lindberg  Jr.  &  Longman, 
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Figure  3.3.  Response  of  SDOF  £  =  0.6 


1984;  Balas,  1978b,  1978a;  Meirovitch  k  Oz,  1978;  Meirovitch  k  Baruh,  1983;  Can- 
field  k  Meirovitch,  1994).  Generally  these  modal  control  methods  belong  either  to 
the  class  of  coupled  methods  or  to  the  class  of  independent  modal  space  control  meth- 
ods developed  by  Meirovitch  and  co-workers  (Meirovitch  &  Oz,  1978;  Meirovitch  k 
Baruh,  1983;  Canfield  k  Meirovitch.  1994).  Using  coupled  methods,  the  closed  loop 
equations  of  the  system  are  coupled  using  feedback  control  wherein  the  optimal  com- 
putation of  the  feedback  gains  require  the  solution  of  a  couple  matrix  Riccati  equation 
(Balas,  1978b,  1978a).  For  large  flexible  structures  the  solution  of  the  Riccati  equation 
can  pose  serious  difficulties  which  limit  significantly  the  applicability  of  the  coupled 
modal  control  methods.  The  IMSC  method,  however,  avoids  such  serious  limitations, 
as  the  control  laws  are  designed  completely  in  the  modal  space,  using  the  uncoupled 
open  loop  equations  of  the  system  as  a  set  of  independent  second  order  equations 
even  after  including  the  feedback  controllers.  Meirovitch  and  co-workers  (Meirovitch 
k  Oz,  1978;  Meirovitch  k  Baruh.  1983;  Canfield  k  Meirovitch,  1994)  showed  that 
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under  such  conditions  it  is  possible  to  compute  the  closed  form  the  optimal  modal 
feedback  gains. 

However,  the  IMSC  method  requires  the  use  of  as  many  controllers  as  the  number 
of  modes  to  be  controlled.  Such  a  requirement  results  in  a  practical  limitation  of  the 
method  when  applied  to  large  structures  where  the  number  of  modeled  modes  can  be 
very  large. 

Lindberg  and  Longman  (1984)  proposed  to  modify  the  IMSC  by  using  a  small 
number  of  controllers  to  control  all  the  modeled  modes  through  the  approximate 
pseudo-inverse(PI)  realization  of  the  modal  controller.  This  modification  can  result 
in  physical  control  forces  which  can  be  far  from  desired  because  the  PI  is,  in  effect,  a 
least  squares  fit  of  N  modal  forces  to  obtain  P  physical  forces.  When  N  equals  P  the 
realization  is  exact  and  is  also  the  same  as  IMSC.  As  N  becomes  much  smaller  than 
P,  the  accuracy  in  obtaining  P  physical  forces  from  N  model  forces  becomes  poor. 
Accordingly,  when  the  realized  forces  are  transformed  back  to  the  modal  space  the 
resulting  modal  forces  will  be  very  far  from  the  optimal  forces  and  this  will  result  in 
deterioration  in  the  performance  of  the  controller. 

For  these  reasons,  the  Modified  Independent  Modal  Space  Control  (MIMSC) 
method  was  initiated  (Baz  &  Poh,  1987).  The  MIMSC  modifies  the  IMSC  algorithm 
to  account  for  the  control  spillover  from  the  controlled  modes  into  the  uncontrolled 
modes  when  a  small  number  of  controllers  is  used  to  control  a  large  number  of  modes. 
This  method  also  incorporates  an  optimal  placement  procedure  for  determining  the 
optimal  location  of  the  controllers  in  the  structure.  Moreover,  the  MIMSC  method 
relies  on  an  efficient  algorithm  for  "time-sharing"  a  small  number  of  controllers  in 
the  modal  space  to  control  a  large  number  of  modes.  The  MIMSC  uses  N  optimally 
placed  controllers  to  control  the  N  modes  that  have  the  highest  modal  energy  at  any 
instant  of  time  and  time  share  these  controllers  among  the  other  residual  modes  when 
the  control  spillover  makes  the  modal  energy  higher  than  the  controlled  modes. 
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3.3,2     Analytical  Results  For  Designing  On-Off  Controller  for  MDOF  systems 

In  this  section,  a  new  design  methodology  is  proposed  for  using  on-off  control 
to  achieve  desired  damping  in  modes  of  interest.  This  method  is  based  on  energy 
transfer  (from  kinetic  to  potential),  as  well  as  on  the  dissipative  nature  of  damping 
forces. 

Time-sharing  of  controllers  was  proposed  in  (Baz  &  Poh,  1987,  1991),  and  these 
ideas  are  used  here. 

Consider  a  second  order  dynamical  system  for  the  n-dimensional  displacement 
vector  x  given  by 

Mx  +  Kx  =  F.  (3.3.12) 

where  M  =  AIT  >  0  is  the  mass  matrix  of  the  structure,  K  =  KT  >  0  is  the 
stiffness  matrix  of  the  structure,  and  F  is  n-dimensional  vector  of  forces  acting  on 
the  structure. 

Equation  (3.3.12)  may  written  in  modal  coordinates  by  using  the  following  weighted 
modal  transformation: 

x  =  $u,  (3.3.13) 

where  u  represents  the  vector  of  modal  coordinates  of  the  system,  and  $  is  the  system's 
modal  shape  matrix  (i.e.  the  matrix  of  eigenvectors). 

Thus,  under  the  transformation  (3.3.13),  equation  (3.3.12)  reduces  to  the  form 

ii  +  Xu  =  f,  (3.3.14) 

where  A  is  the  diagonal  matrix  of  eigenvalues  of  M~'  K  matrix  and  /  =  $rF. 
The  controller  being  on,  yields  the  dynamics: 

Mx  +  Kx  =  -AA'i.  (3.3.15) 
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where  (K  +  AA)7'  =  (K  +  AA)  >  0. 

Under  switching  control,  it  is  desired  that  the  system  behave  like  a  damped  dy- 
namical system  of  the  form 

Mx  +  Di  +  Kx  =  0,  (3.3.16) 

where  D  is  an  n  x  n  desired  proportional  damping  matrix:  i.e.  D  =  aM  +  iJK,  where 
a,0  are  scalar  constants.  Note  that  the  mass  and  stiffness  matrices  are  the  same. 

The  energy  stored  (or  released)  in  time  T  due  to  the  added  stiffness  at  the  con- 
troller is  given  by 

E±K=l-fa    jt{xT\Kx)dt.  (3.3.17) 

The  energy  dissipated  due  to  the  damping  term  D  over  an  interval  of  length  T  is 
given  by 

E^\lTJt(xTm)dt-  (3-3'18) 

Assuming  the  controller  is  only  on  for  T  seconds,  and  using  equations  (3.3.13)  and 
(3.3.17)  the  contribution  of  the  ith  mode  to  the  energy  stored  in  the  controller  is 

E±Ki  =  \[  jf  (u,  [$f  AK$i]  u.)  dt.  (3.3.19) 

where  $,  is  the  ilk  column  of  5>.  The  contribution  of  the  ith  mode  to  the  energy  which 
would  be  dissipated  in  the  desired  dynamics  through  the  damping  over  the  period  T; 

is 

ED,  =  \f'jt  («i  [*f£*i]  «i)  dt.  (3.3.20) 

Here  T,  is  taken  as  the  period  of  the  iih  mode  of  vibration. 
Since  proportional  damping  is  assumed. 

*f  £>*;  =  2QuJi,  (3.3.21) 
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where  Q  is  the  desired  modal  damping  and  Wj  is  the  natural  frequency  of  the  itk  mode. 
If  Ui  is  assumed  to  be  of  the  form  u,  =  e'^'sin  fat),  then 

u((t)  =  -Qe~(itsm  fat)  +  e'^facos  fat)  (3.3.22) 

and 

&i(t)  =  C,2e"c,'sm  fat)  -  2Qe~i'facos  fat)  -  e"c<*w2s»n  fat) .  (3.3.23) 

The  structure  does  work  on  the  controller  during  only  half  of  the  cycle  of  vibration 

of  a  given  mode  and  during  the  other  half  of  the  cycle,  the  controller  releases  energy 

into  the  structure.   This  logic  motivates  the  scheme  for  switching  the  controller  on 

only  for  half  of  the  period  of  the  mode  to  be  controlled.    Using  T  =  \  the  two 

integrals  (in  equations  (3.3.18)  and  (3.3.19))  may  be  evaluated,  and  the  contribution 

of  mode  i  to  the  energy  stored  in  the  controller  over  its  period  (when  the  controller 

is  only  on  for  half  the  period)  is  equated  to  the  contribution  of  mode  i  to  the  energy 

dissipated  in  the  damping  term  over  its  period.   This  yields  an  expression  for  AK 

(nonzero  elements  of  AK).   Although  there  is  no  restriction  that  AK  be  diagonal, 

for  simplifying  the  expressions  (and  to  make  them  more  clear),  AK  is  assumed  to  be 

diagonal. 

Thus 

$£A#$m=  Y.  Qjm&Kjfijn.  (3.3.24) 

j  such  that  &Kjj^0 

This  depends  only  on  one  those  entries  in  the  mlh  column  of  $,  for  which  the  corre- 
sponding entries  on  the  diagonal  matrix  AK  are  non  zero. 

Take  first  the  case  that  only  one  entry  in  AK  is  nonzero.  Let  j  mark  the  index 
of  this  entry. 

Evaluating  the  derivatives  in  the  integrals  of  equations  (3.3.19)  and  (3.3.20)  (eval- 
uating the  first  integral  for  only  half  the  cycle  as  described  earlier)  gives 

i  rTi/2 


Eak'  =  J /       2  (^^AATjjtti)  dt  (3.3.25) 
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and 

1    fli 


ED.  =  -  I '  '  fuj2CiW,Ui  +  u?2C,«*1  dt.  (3.3.26) 

2  Jo     y  ' 

Using  the  expressions  given  earlier  for  u,  it.  u,  $f  AA'$;,  $J D<b(,  given  by  equa- 
tions (3.3.22.3.3.23.3.3.24,3.3.21)  define  gl;,.^,.^,  (the  expanded  form  of  the  inte- 
grand in  equation  (3.3.25)  and  terms  in  the  integrand  of  equation  (  3.3.26))  as 

glt  =  [""  (-C,e_Ci'sm(wjt)  +  e'':'t^1cos(ujtt))e-c'tsin(uitt)AKjj^2j)dt        (3.3.27) 


cj2,  =  /   '  2(C,?  e-'-itsin(uJit)+2(-C,i)e-^tuiicos(ujtt)-e-iitu)'2isin(ujit))e''Qtsin(ujit))  Q  ujtdt. 
Jo 

(3.3.28) 
53;=/   '  2(-Cie~c''sm(o;ii)  +  e-c"'u,- cos(wii))2u;,Cid*-  (3.3.29) 

JO 

Since  we  are  considering  the  contribution  of  only  one  mode  at  this  point,  there 
will  be  a  nonzero  contribution  from  only  one  diagonal  element  of  AK  matrix.  Solving 
gU  =  g2i+g3i  gives  an  expression  for  the  nonzero  diagonal  element  in  the  AK  matrix 
as 

AKn  =  -i!±£  (3.3.30) 

u 

When  there  is  more  than  one  mode  to  control,  the  matrix  AA'  should  have  only 

as  many  non-zero  diagonal  entries  as  the  number  of  modes  to  be  controlled.  In  this 

case,  gli  takes  the  form: 

J2  f"'  (-Qe-Qtsin{uj,t)  +  e~(-'tuicos(LOit))e<'tsm(uiit)AK,j^2ijdt. 

]  such  that  ilA'jjjtIK0 

In  such  cases,  the  integrand  in  equation  (3.3.25)  and  equation  (3.3.26)  will  have 
as  many  terms  as  the  number  of  modes  to  be  controlled  and  as  many  equations  of 
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the  form  3I,  =  </2;  +  fl3,  as  are  necessary  will  have  to  be  solved  in  order  to  obtain 
values  of  different  non-zero  elements  of  the  AA"  matrix.  With  more  than  one  mode 
to  control,  there  will  be  m  expressions  of  the  form  (3.3.25)  and  m  expressions  of  the 
form  (3.3.26)  (one  for  each  of  the  desired  modes  to  be  controlled).  This  will  lead 
to  m  equations  of  the  form  (3.3.27)-(3.3.29)  which  can  be  solved  for  the  m  nonzero 
diagonal  entries  of  AA'  matrix. 

The  flowchart  for  the  proposed  MDOF  on-off  control  scheme  is  given  in  Figure 
3.4.  Note  that  at  any  given  time,  controllers  are  only  set  to  act  on  one  mode  in  order 
to  increase  its  damping.  The  switching  procedure  computes  the  energy  in  different 
modes  of  the  structure  and  checks  whether  the  energy  in  one  of  the  desired  modes  is 
the  highest.  If  it  finds  that  energy  in  one  of  the  desired  modes  is  taking  the  highest 
value,  for  half  the  period  of  vibration  of  this  mode,  all  the  controllers  are  dedicated 
to  damp  the  vibrations  in  this  mode  (the  timing  is  kept  track  by  the  use  of  a  counter 
Tc).  During  this  half  cycle  of  vibration,  controllers  are  individually  switched  on  and 
off  depending  upon  whether  or  not  work  is  done  on  them  by  the  structure. 

3.3.3     Four  Degree  of  Freedom  example  with  damping  in  one  mode 

Consider  a  four  degree  of  freedom  undamped  system  with  mass  and  stiffness  prop- 
erties 


Mass  = 


2 

0 

0 

0  1 

0 

3 

0 

0 

0 

0 

4 

0 

0 

0 

0 

5 

(3.3.31) 


Stiffness  = 


3-200 
-2      4-2      0 
0-2      5-2 
0       0-26 


(3.3.32) 
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The  controller  was  placed  at  the  first  mass  and  the  damping  in  the  first  mode  was 
desired  to  be  0.2.  Unlike  the  third  example  (Kabe's  problem)  in  this  section,  this 
model  had  mass  normalized  eigenvector  components  of  the  same  order  (components 
of  different  eigenvectors  at  nodal  locations).  Equation  3.3.30  suggests  that  controller 
gains  are  inversely  proportional  to  the  eigenvector  components.  This  suggested  that 
placement  of  controller  at  any  of  the  masses  would  have  resulted  in  a  value  of  AA' 
in  the  same  order.  The  on-off  controller  scheme  proposed  in  Section  4.3.2  was  im- 
plemented. Nodal  and  Modal  time  responses  of  the  system  are  given  in  Figures  3.5 
and  3.6  respectively.  From  the  plot  of  modal  responses,  it  is  clear  that  the  control 
scheme  introduces  damping  in  first  mode.  The  value  of  this  damping  is  found  to  be 
0.19.  Figure  3.7  gives  the  energy  in  different  modes  and  Figure  3.8  gives  the  velocity 
and  modal  force  for  the  first  mode  of  vibration.  From  Figure  3.8  it  is  clear  that  the 
control  scheme  introduces  damping  into  first  mode. 

3.3.4     Four  Degree  of  Freedom  example  with  damping  in  all  modes 

In  this  example,  controller  parameters  used  in  the  simulation  are  such  that  all  four 
modes  of  the  four  DOF  system  given  in  equations  (3.3.31)  and  (3.3.32)  exhibit  the 
desired  damping.  The  desired  damping  values  in  the  four  modes  are  .2, .15.0.20  and 
0.17.  The  natural  frequencies  of  the  uncontrolled  system  are  0.556.  0.9657,  1.2791  and 
1.5511  radians.  The  system  was  simulated  and  the  algorithm  given  in  the  previous 
section  was  implemented.  It  was  observed  that  system  behaves  like  a  well  damped 
system  with  modal  damping  very  close  to  the  desired  values.  This  was  simulated  with 
several  initial  conditions.  Figures  3.9-3.12  give  the  nodal  response  and  Figures  3.13- 
3.16  give  the  modal  response  of  this  system.  The  observed  modal  damping  values  are 
0.18,  0.14,  0.18  and  0.16.  The  natural  frequencies  of  the  system  with  control  acting 
upon  it  are  identified  as  0.7454,  0.9827,  1.1310  and  1.2454  radians/second  respectively, 
matching  exactly  with  the  natural  frequencies  of  the  uncontrolled  system. 
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3.3.5     Kabe's  Problem  -  With  damping  in  only  two  modes 

Kabe's  eight  degree  of  freedom  model  is  shown  in  Figure  3.17.  The  natural  fre- 
quencies of  this  uncontrolled  system  are  30.0280,  31.7163,  31.7557.  32.3020,  33.3267, 
35.5763.  38.7887  and  41.8884  radians/second  respectively.  The  mass-normalized 
eigenvectors  of  this  system  are  listed  in  Table  4.1.  In  order  to  avoid  high  values 
of  controller  gains,  the  controllers  were  placed  at  4"1  and  6"1  masses.  The  controller 
locations  were  selected  based  upon  the  mass  normalized  eigenvector  components.  The 
proposed  design  procedure  shows  that  the  reciprocal  of  the  square  of  the  eigenvector 
component  at  the  location  of  the  controller  decides  the  gain  gain  of  the  controller  (i.e. 
if  damping  is  desired  in  mode  j  and  controller  is  placed  at  location  g,  reciprocal  of 
the  square  of  the  gih  component  of  the  jth  mass  normalized  eigenvector  decides  the 
controller  gain).  Note  that  the  fourth  and  sixth  mass  normalized  eigenvector  have 
fourth  and  sixth  components  which  are  within  reasonable  limits  (leading  to  reason- 
able values  for  the  elements  of  AK  matrix.  The  desired  damping  in  the  fourth  and 
sixth  mode  were  set  to  0.10  and  0.15,  respectively.  The  initial  conditions  were  set 
such  that  the  initial  energy  was  concentrated  in  the  fourth  and  sixth  modes.  When 
these  conditions  were  simulated,  it  was  seen  that  the  control  algorithm  achieves  the 
desired  goals  within  reasonable  limits.  Figures  3.18-3.25  give  the  nodal  response  and 
Figures  3.26-3.25  give  the  modal  response  of  this  system.  In  order  to  show  the  impor- 
tance of  the  controller  locations,  the  same  system  was  simulated  with  the  controllers 
located  at  the  first  and  the  second  mass.  The  desired  damping  in  fourth  mode  and 
sixth  mode  were  set  to  0.1  and  0.15.  Controller  gains  in  this  case  were  42  and  38166 
compared  to  2724  and  3486  in  the  earlier  case.  Dampings  achieved  in  this  case  were 
.06007  and  0.1076  in  mode  4  and  6  compared  to  the  values  of  0.08  and  0.14  in  the 
previous  case.  In  order  to  implement  the  proposed  scheme,  the  controller  gains  would 
have  to  be  such  that  they  do  not  significantly  affect  the  mass  of  the  system.   If  the 
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size  and  weight  of  the  controllers  change  the  mass  of  the  system  the  effects  of  this 
will  have  to  be  taken  into  consideration. 


Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Mode  6 

Mode  7 

Mode  8 

-2.1170 

-30.0002 

-2.5940 

-10.6279 

-5.6896 

-2.6102 

-31.6226 

0.0002 

-0.0267 

-0.3125 

-0.0269 

-0.1023 

-0.0467 

-0.0129 

0.0030 

0.00001 

-0.2832 

-0.0319 

0.0040 

0.3769 

0.4901 

0.3401 

-0.0001 

0.0001 

-0.5873 

-0.0020 

0.0068 

0.2612 

0.0014 

-0.5282 

0.0001 

-0.0005 

-0.5973 

0.0329 

-0.0806 

-0.6016 

-0.0013 

0.1988 

0.0001 

0.0001 

-0.2841 

-0.0029 

0.0827 

0.3724 

-0.4890 

0.3361 

0.0001 

0.0032 

-0.0308 

-0.0085 

0.6378 

-0.1347 

0.0598 

-0.0172 

0.0001 

0.0020 

-8.0992 

-0.2782 

16.9175 

8.5882 

-15.5332 

14.9251 

-0.0005 

-22.3605 

Table  3.1.  Mass  normalized  eigenvectors  of  Kabe's  model 
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INPUT  STRUCTURAL  PARAMETERS,  CONTROLLER  LOCATION^ 
AND  LOADING 


Set  Controller  flap 
off.  Reset  To- 


J 


Compute  time  history  of  all 
nodes  in  physical  coordinates 


Compute  time  history  of  all 
nodes  in  modal  coordinates 


Switch  the  controllers  on  if  the 
work  done  on  them  is  positive. 


Compute  energy  in  all  modes. 
Find  max  Energy  in  desired  mode  j 


Store  j. 

Set  Controller  flag  on. 


Figure  3.4.  Flowchart  of  MDOF  damping  control  algorithm 
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Figure  3.5.  Nodal  Response  of  a  4-DOF  System  with  (ides  =  0.2 
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Figure  3.6.  Modal  Response  of  a  4-DOF  system  with  Cides  =  0.2 
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Figure  3.7.  Energy  in  different  modes  of  a  4-DOF  system  with  (lrfes  =  0.2 
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Figure  3.8.  Modal  force  and  velocity  (for  the  first  mode) 
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Figure  3.9.  Response  at  Node  1  of  a  4-DOF  system  (with  all  modes  controlled) 


Figure  3.10.  Response  at  Node  2  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.11.  Response  at  Node  3  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.12.  Response  at  Node  4  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.13.  Mode  1  Response  of  a  4-DOF  system  (with  all  modes  controlled) 


Figure  3.14.  Mode  2  Response  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.15.  Mode  3  Response  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.16.  Mode  4  Response  of  a  4-DOF  system  (with  all  modes  controlled) 
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Figure  3.17.  Kabe's  Model 
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Figure  3.18.  Controlled  response  at  Node  1  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.19.  Controlled  response  at  Node  2  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.20.  Controlled  response  at  Node  3  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.21.   Controlled  response  at  Node  4  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.22.  Controlled  response  at  Node  5  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.23.  Controlled  response  at  Node  6  of  Kabe's  Model  (with  only  two  modes 
controlled) 


52 


lime  in  seconds 


Figure  3.24.  Controlled  response  at  Node  7  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.25.  Controlled  response  at  Node  8  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.26.    Controlled  mode  1  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.27.    Controlled  mode  2  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.28.    Controlled  mode  3  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.29.    Controlled  mode  4  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.30.    Controlled  mode  5  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.31.    Controlled  mode  6  Response  of  Kabe's  Model  (with  only  two  modes 

controlled) 


time  in  seconds 


Figure  3.32.    Controlled  mode  7  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 
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Figure  3.33.    Controlled  mode  8  Response  of  Kabe's  Model  (with  only  two  modes 
controlled) 


3.3.6     Kabe's  Problem  -  With  damping  in  all  modes 


In  this  section  Kabe's  eight  degree  of  freedom  system  was  simulated  with  damp- 
ing introduced  in  all  modes.  The  desired  damping  values  in  the  eight  modes  are 
.20,0.18,0.15.0.10.0.10.0.15.0.13  and  0.12.  Figures  3.34-3.41  show  the  modal  responses 
of  this  system.  From  these  figures  it  is  clear  that  controlling  all  the  modes  of  the 
system  gives  better  performance  when  compared  co  the  case  with  only  few  modes  con- 
trolled. The  achieved  damping  values  in  the  eight  modes  are  0.19,0.16.0.14.0.08,0.14.0.14,0.11 
and  0.11.  The  oscillations  that  are  seen  in  Figure  3.41  are  due  to  the  activation  of  the 
local  modes  of  vibration.  The  natural  frequencies  of  the  controlled  modes  were  iden- 
tified as  30.0280.  31.7163,  31.7557.  32.3020,  33.3267.  35.5763.  38.7887  and  41.; 
radians/second  respectively,  which  match  the  uncontrolled  natural  frequencies. 
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Figure  3.34.  Controlled  mode  1  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.35.  Controlled  mode  2  Response  of  Kabe's  Model(with  all  modes  controlled) 
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Figure  3.36.  Controlled  mode  3  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.37.  Controlled  mode  4  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.38.  Controlled  mode  5  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.39.  Controlled  mode  6  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.40.  Controlled  mode  7  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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Figure  3.41.  Controlled  mode  8  Response  of  Kabe's  Model  (with  all  modes  controlled) 
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3.4     Summary 

A  design  methodology  for  the  energy  based  on-off  switching  decentralized  control 
of  mechanical  structures  is  presented.  Energy  consumption  is  minimized  by  turning 
on  the  controllers  when  a  structure  is  doing  work  on  the  controllers,  and  turning  them 
off  after  prescribed  on-times.  For  a  SDOF  system,  closed  form  solution  is  given  for 
designing  the  controller  parameter.  For  a  MDOF  system  the  proposed  methodology 
can  be  used  to  achieve  a  desired  average  damping  rate  when  each  mode  may  be 
uniquely  controlled.  The  controller  design  methodology  and  the  switching  schemes 
are  presented.  Though  the  proposed  scheme  suggests  the  use  of  no  external  energy, 
in  practice,  such  methods  would  require  some  energy  input  for  signal  monitoring,  as 
well  as  to  drive  an  electro-mechanically  switched  controller  using  a  computer  or  a 
microprocessor. 

Since  altering  the  configuration  of  sensors  and  controllers  does  not  alter  the  fact 
that  structures  with  positive  damping  dissipate  energy,  the  passivity  theorem  (Morris 
&  Juang,  1994)  can  be  applied  in  this  case  to  prove  that  modal  forces  from  the  con- 
trollers will  not  destabilize  the  system.  Knowledge  of  the  system  is  used  in  obtaining 
the  controller  parameters.  Changes  or  errors  in  the  system  or  controller  parameters, 
however,  will  not  destabilize  the  system  because  it  is  assumed  to  be  dissipating  energy. 


CHAPTER  4 
FUEL-EFFICIENT  WAYS  OF  VIBRATION  CONTROL 

4.1     Introduction 

The  work  on  switching  control  presented  in  Chapter  2  and  3  was  motivated  by  a 
study  of  the  literature  which  used  on-off  control  of  vibrations,  such  as  Athanassiades 
(1963).  Athans  (1964b,  1964a),  Foster  and  Silverberg  (1991),  Silverberg  (1986),  Red- 
mond and  Silverberg  (1992)  and  Medith  (1964).  The  work  of  Foster  and  Silverberg 
(1991)  is  one  such  study  which  uses  fuel  minimization  as  an  optimization  criterion  for 
on-off  control  of  structural  vibrations.  In  that  work,  the  minimum  fuel  use  control 
law  imposed  that  control  pulses  be  applied  every  time  a  zero  crossing  of  the  the  po- 
sition (and  maximum  value  of  the  velocity)  would  occur.  A  plot  of  fuel  consumption 
versus  sampling  time  was  also  presented  in  (Foster  &  Silverberg,  1991).  This  chapter 
is  devoted  to  the  study  of  the  effects  of  sampling  time  and  window  size  on  the  fuel 
consumption. 

4.2     Minimum  Fuel  Control  of  Mechanical  Systems 

Foster  and  Silverberg  (1991)  developed  a  minimum  fuel  based  switching  control 
law  for  MDOF  second  order  systems  of  the  form: 

Mx  +  Kx-F,  (4.2.1) 

where  M  —  M'  is  the  n  x  n  mass  matrix  of  the  structure,  K  is  the  n  x  n  stiffness 
matrix  of  the  structure,  x  and  x  are,  respectively,  the  n-dimensional  displacements 
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and  accelerations  at  the  nodal  points  of  the  structure,  and  F  is  the  n-dimensional 
vector  of  forces  acting  on  the  structure. 

Since  the  control  law  depends  on  points  (in  time)  of  maximum  velocity  (or  mini- 
mum position)  these  are  determined  using  running  averages  of  velocities  and  positions: 

fif(t)  =  J  a{s,t)f(s)ds,    J  a(s,t)ds  =  1  (4.2.2) 

ij(t)  =  £  a(s,t)U(s)  -  HfWfds,  (4.2.3) 

where  a(s,  i)  is  a  window  function  of  the  form: 


,    ,     .  l/r,  s>t-T  , .  n .. 

a(s,t)  =  {     '0i       g<t_T  (4.2.4) 


Here,  T  is  the  window  size. 

Equations  (4.2.2)  and  (4.2.3)  are  used  in  determining  the  control  weightings  for  a 
piecewise  constant  (vector)  control  law. 

In  that  work,  the  fuel  consumption  C(t)  in  time  t  is  measured  using 

Hi)  =  E   /'l/.l  ds  (4.2.5) 

where  /4(s)  is  the  ith  component  of  the  control  force  at  time  s,  in  modal  coordinates. 

4.3     Effects  of  changes  in  window  size  and  sampling  period 
4.3.1     Example  1 

Foster  and  Silverberg  (1991)  worked  out  an  example  prescribed  by  the  following 
parameters: 

M  =  1,  K  =  1  and  T  =  10. 
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where  T  is  the  number  of  samples  in  the  window.  In  Figure  4.1  the  fuel  consump- 
tion was  plotted  as  a  function  of  number  of  samples  per  period  of  oscillation  (denoted 
by  N,  a  parameter  which  determines  the  sampling  period).  Notice  the  sharp  increase 
in  fuel  consumption  for  certain  values  of  N. 

As  it  turns  out,  when  these  increases  in  fuel  consumption  occur,  so  does  a  trig- 
gering of  the  higher  harmonics  in  the  control  pulses.  The  fast  fourier  analysis  of  the 
control  pulses  shows  that  the  ratio  of  signal  power  at  the  fundamental  frequency  to 
the  sum  of  the  signal  power  at  the  higher  harmonic  frequencies  decreases  at  the  points 
of  positive  slope  in  this  plot. 


20  30  40 

Number  of  samples  per  period  of  oscillation 

Figure  4.1.  Effect  of  Number  of  Samples  per  period  of  oscillation  on  fuel  consumption 
(example  1) 


The  frequency  spectra  of  the  power  of  control  signal  are  shown  in  Figures  4.2-4.5 
for  selected  points  on  the  plot  shown  in  Figure  4.1.  The  fundamental  frequency  of  this 
system  is  0.159  Hz.  The  ratio  of  the  power  in  the  control  signal  at  the  fundamental 
frequency  to  that  represented  by  the  higher  harmonics  decreases  at  the  points  of 
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positive  slope  on  the  plot  in  Figure  4.1.  This  is  seen  in  the  plots  in  4.3  and  4.5.  In 
these  plots  note  that  for  values  of  N  =  15  and  N  =  25,  the  power  at  higher  harmonics 
take  higher  values  when  compared  to  plots  for  values  of  N  where  the  slope  of  the  fuel 
curve  shown  in  Figure  4.1  is  negative.  At  the  points  of  positive  slope  of  the  plot  in 
Figure  4.1.  sampling  time  is  such  that  the  window  size  becomes  an  integral  multiple 
of  one  of  the  higher  harmonics.  Note  that  for  N  =  15  window  size  is  twice  the  period 
of  third  harmonic  and  for  N  =  25  window  size  is  twice  the  period  of  fifth  harmonic. 


0.2  0.3  0.4  0.5  0.6  0.7 

Frequency  in  Hz  (Number  ot  samples=1 1 ) 


Figure  4.2.  Frequency  spectrum  of  control  pulse  in  example  1  (N=ll) 
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0.4  0.6  0.8 

Frequency  in  Hz  (Number  of  samples=15) 


Figure  4.3.  Frequency  spectrum  of  control  pulse  in  example  1  (N=15) 
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0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8 

Frequency  in  Hz  (Number  ol  samples=21 ) 

Figure  4.4.  Frequency  spectrum  of  control  pulse  in  example  1  (N=21) 

4.3.2     Example  2 

In  this  example,  the  system  parameters  used  are:  M  =  1,  K  =  2,  and  T  =  10. 

The  system  was  simulated  using  the  same  control  scheme  and  the  effects  of  sam- 
pling time  on  fuel  consumption  was  studied.  The  plots  of  fuel  consumption  and 
spectra  of  control  signal  are  given  in  Figures  4.6-4.10.  The  natural  frequency  of  vi- 
bration is  .23  Hz.  Just  like  in  the  previous  example,  the  higher  harmonics  in  the 
control  signal  make  the  fuel  consumption  to  increase  at  certain  values  of  sampling 
times. 

4.4     Summary 


The  effects  of  sampling  on  fuel  consumption  using  Foster  and  Silverberg  (1991)'s 
minimum  fuel  on-off  control  law  was  studied.  It  was  noted  that  when  the  significant 
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0.016 


0.2  0.4  0.6  0.8  1  1.2  1.4  1.6  1.8  2 

Frequency  in  Hz  (Number  of  samples=25) 

Figure  4.5.  Frequency  spectrum  of  control  pulse  in  example  1  (N=25) 


harmonics  in  the  control  signal  impose  that  window  size  is  integer  multiple  of  the 
period  of  higher  harmonics,  fuel  consumption  was  increased. 


Tl 


10  20  30  40  50 

Number  of  samples  per  period  of  oscillation 
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Figure  4.6.  Effect  of  Number  of  Samples  per  period  of  oscillation  on  fuel  consumption 
(example  2) 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.E 

Frequency  in  Hz  (Number  of  samples=9) 

Figure  4.7.  Frequency  spectrum  of  control  pulse  in  example  2  (N=9) 
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0.2  0.4  0.6  0.8 

Frequency  in  Hz  (Number  of  samples=11) 


Figure  4.8.  Frequency  spectrum  of  control  pulse  in  example  2  (N=ll) 


I  0.008 


0.5  1  1.5  2  2.5 

Frequency  in  Hz  (Number  of  samples=25) 

Figure  4.9.  Frequency  spectrum  of  control  pulse  in  example  2  (N=25) 
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1  1.5  2  2.5  3  3.5 

Frequency  in  Hz  (Number  of  samples=38) 


Figure  4.10.  Frequency  spectrum  of  control  pulse  in  example  2  (N=38) 


CHAPTER  5 
CONCLUSION  AND  SUGGESTIONS  FOR  FUTURE  WORK 

5.1     Discussion 

This  study  proposed  new  schemes  for  structural  control  design  using  the  energy 
input  into  the  controller  to  damp  out  the  vibrations.  Through  simulations  it  was 
shown  that  such  a  decentralized  control  law  can  be  efficient  in  improving  the  perfor- 
mance of  structures.  For  SDOF  systems,  the  proposed  design  methodology  introduces 
desired  damping  within  very  close  bounds.  For  MDOF  systems,  using  the  techniques 
of  Independent  Modal  Space  Control,  a  design  algorithm  was  developed  to  introduce 
specified  damping  into  desired  modes  of  vibration.  The  controllers  were  decentralized 
in  that  the  decision  for  switching  a  given  controller  was  based  on  feedback  from  local 
states  only.  It  was  also  shown  how  the  uncontrolled  modes  are  affected  by  controlling 
only  few  modes  in  multi  degree  of  freedom  systems.  This  suggests  that  the  proposed 
methodology  be  used  to  control  all  of  the  dominant  modes  of  the  vibrating  structure. 

5.2     Future  Research 

The  minimum  energy  on-off  control  methods  proposed  in  this  work  should  increase 
the  life  span  of  space  structures  since  very  minimal  energy  is  needed  from  the  external 
power  with  such  a  method. 

To  fully  demonstrate  the  practical  advantage  of  these  schemes,  one  should  de- 
termine the  kind  of  controllers  that  can  be  used  for  such  applications.     A  simple 
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gyroscope,  magnetic  bearing  gyros,  or  torque  wheels  provide  potentially  viable  op- 
tions for  such  applications.  Control  momentum  gyros  have  been  effectively  used  for 
the  attitude  control  of  spacecrafts.  For  instance,  single  gimbaled  gyros  have  been 
installed  on  the  Soviet  MIR  station,  while  double  gimbaled  gyros  were  used  in  the 
NASA  Skylab.  Momentum  gyro-based  control  designs  have  been  tried  in  the  recent 
past  for  slewing  maneuvers  and  vibration  control  of  space  structures.  When  used  in 
conjunction  with  the  proposed  on-off  control  scheme,  momentum  gyros  might  pose 
additional  problems  associated  with  the  controller  dynamics  due  to  the  delays  in- 
troduced bv  them.  The  proposed  methodology  is  based  on  the  assumption  that  the 
controllers  that  are  used  do  not  release  the  absorbed  energy  back  into  the  structure. 
The  energy  absorbed  by  the  controller  can  be  either  stored  or  dissipated  into  some 
sink  (like  a  heating  element).  Peizo-electric  crystal  based  controllers  can  convert  the 
force  into  electricity  which  may  be  stored.  Gyro  based  controllers  can  convert  the 
forces  into  angular  momentum  which  may  be  used  to  drive  small  generators  to  gener- 
ate electrical  energy.  The  effects  of  such  controllers  on  the  mass  of  the  structure  will 
have  to  be  investigated. 
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